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The non- normality of Wilson-type lattice Dirac operators has important consequences — the application of 
the usual concepts from the textbook (hermitian) quantum mechanics should be reconsidered. This includes an 
appropriate definition of observables and the refinement of computational tools. We show that the truncated 
singular value expansion is the optimal approximation to the inverse operator D -1 and we prove that due to 
the 75-hermiticity it is equivalent to 75 times the truncated eigenmode expansion of the hermitian Wilson-Dirac 
operator. 



1. INTRODUCTION 

Wilson-type lattice Dirac operators are non- 
normal, i.e. 

[D.D^^O, (1) 

and this fact has several important consequences. 
The essential difference to normal operators is 
that they are not diagonalizable by a unitary 
transformation. Instead they are diagonalizable 
(if at all) by a similarity transformation of the 
form 

diag(A 1 ,...,A Ar ) = X- 1 DX, (2) 

with N = dim(Z3). It follows that for a given 
eigenvalue there are distinct left (Li | and right 
\Ri) eigenvectors 

D\IU) = Xi\Ri) , (3) 

(L i \D = X i (L i \, (4) 

in the sense that (in general) (L,| ^ (|i?,;))^. 

This fact must be taken into account when per- 
forming the spectral decomposition 

JV JV 

D = J2*i\Ri)(Li\¥=Y, X i\ R iK R i\' ( 5 ) 
i=l i=l 

but also raises the question of an appropriate gen- 
eralization of the other concepts familiar from the 
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textbook (hermitian) quantum mechanics. For 
example, if an observable in the continuum is de- 
fined as (ipi\A\ipi) , where \ipi) is an eigenstate of 
the continuum operator D — what is the proper 
generalization for the case where the left and the 
right eigenvectors are distinct? 

In we advocate the use of (Li\A\Ri) instead 
of (Ri\A\Ri) as a kind of improved observable. 
Note that in lattice field theory this question is 
not of importance in the continuum limit, because 
the difference between the two definitions disap- 
pears. However, this question was also raised in 
the domain of solid state physics — see Q for an 
interesting discussion. 

2. TRUNCATED EIGENMODE EXPAN- 
SION 

In this contribution we want to examine a dif- 
ferent problem, which is also connected to the 
non-normality of the Wilson-type lattice Dirac 
operators. The idea to approximate the inverse 
operator Z? _1 by just a few eigenmodes with 
eigenvalues close to zero is physically appeal- 
ing. One hopes that the lowest modes dominate 
physics and that the truncated eigenmode expan- 
sion 

fc«JV 

D k 1= E ttI^X^I ( 6 ) 

i—1 
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Figure 1. Correlation function for several values 
of At (At increases as one steps down from the 
top curve) plotted versus the number of eigen- 
modes used to approximate the inverse of the 
Wilson-Dirac operator — for a given k, the cor- 
relation function is computed from , as de- 
fined in (||). The dimension of D is 512, so the 
final points give the values for the exact inverse 
of D -1 . 



with the ordering 
11 1 1 

will be a fair approximation of the proper inverse 

However, in p] it was found that such an ex- 
pansion is unstable and that it does not uniformly 
saturate by increasing k. This is illustrated in 
Fig. [l] where the values of the correlation func- 
tion for different temporal distances At are plot- 
ted versus k. (The illustration is just for two 
dimensional QED on a 16 2 lattice for (3 = 2. 
but it should be stressed that for the discussion 
which follows this is irrelevant — the arguments 
are pure linear algebra and they will be valid for 
any physical system and parameter set.) 

By a closer examination of pseudospectra (for 
more details about pseudospectra see |J) of D, 
we identified that the spikes in Fig. [I] correspond 
essentially to the instabilities due to the non- 
normality of D. For non-normal operators the 
truncated eigenmode expansion in the original 



biorthogonal basis seems to be problematic. How- 
ever, one can do better than that. 

3. SINGULAR VALUE DECOMPOSI- 
TION 

For any matrix A there exists a so-called sin- 
gular value decomposition (SVD) 



A=UY,V ] 



Pi\Ui){Vi 



(8) 



where U and V are unitary and £ is the diagonal 
matrix with positive or zero entries on the diag- 
onal — the "singular values" <7j. We assume the 
following ordering: 



a i — °2 — ■ ■ ■ 17 n — • 



(9) 



Among other interesting properties, one can 
prove the following theorem (Th. 2.5.3. in Q): 

If k < r = rank(.A) and 



k 

E< 



A k = y ]ai\ui)(vi 



(10) 



then 



min \\A-B\\ 2 = \\A-A k \\ 2 = a k+1 , (U) 

rank(B) — k 

i.e., Ak is the closest matrix to A that has rank 
k. In other words, the truncation of SVD is the 
optimal approximation to the original matrix for 
a given k. 

Actually, the drawback is that for a general ma- 
trix A one needs two sets of vectors: |u)-s and \v)- 
s. However, the situation is better for the Wilson- 
type lattice Dirac operators (and their inverses). 
Using the 75-hermiticity property 

= 75^75 (12) 

we are going to prove that the truncated SVD of 
D^ 1 is, identical to 



A 1 



i=l 



(13) 



with the ordering according to the absolute values 
of inverse fi-s: 



1 1 

M ~ IM2I ~ 
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(14) 
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Figure 2. Correlation function for several tempo- 
ral distances At plotted versus the truncated sin- 
gular value expansion, as defined by (|ll|). This 
time, the number of eigenmodes refers to the 
modes of the hermitian Wilson-Dirac operator H 
used to construct H7 1 . 



The fi-s and \w)-s are the eigenvalues and eigen- 
vectors of the hermitian Wilson-Dirac operator 



N 



H = 75-D = WMW^ = ^2 / fjL i \w i ){w. 

l 



(15) 



The hermiticity of H follows from the 75- 



hermiticity property (12). Thus H is diagonaliz- 
able by a unitary transformation W and has real 
eigenvalues /i^ , ordered on the diagonal of the di- 
agonal matrix M. 

Per construction, D'D is a hermitian matrix 
and it is easy to see that 

D ] D = WMW^WMW^ = WM 2 Wl (16) 

However, by using the singular value decomposi- 
tion (D) of D it follows that 

D^D = VYXJ^UYV^ = VT, 2 Vl (17) 

As V is unitary and E 2 diagonal, one must have 
(with appropriate ordering) 

(18) 



Mi 



i.e., the singular values of D are just the absolute 
values of the eigenvalues of the hermitian Wilson- 
Dirac operator H . Also, the |u)-s should corre- 
spond to the \w)-s (there is an arbitrary phase 



factor, but it cancels in (|13|)). Finally, we can 
rewrite: 

D = UT,V^ = Usign{M)MW^ = lb WMW^ (19) 

which makes it obvious that 

U = f 5 Wsiga(M) \ Ui ) =sigmjU i )75k)-(20) 

Hence, the truncated singular value expansion of 
D is 

k k 

D k = / J Vi\ui){vi\ = ^ en sign(/^ ) 75 1 w l ) ( Wj \ 



i=l 



i=l 



: 75 



y~lHi\wi){wi\ = -f 5 Hk ■ 



(21) 



As D 1 also satisfies the 75-hermiticity property 
(p^), the same reasoning applies and (|l^) follows. 

4. CONCLUSION 

Fig. ^illustrates that, as expected, by using the 
truncated singular value expansion, non-normal 
artefacts disappear and the saturation as a func- 
tion of k is superior to the expansion in the eigen- 
modes of D. This approach has already been suc- 
cessfully used in a realistic QCD setting || . 
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